One- and two-neutron removal reactions from the most neutron-rich carbon isotopes 
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The structure of 



^^C has been investigated using high-energy (about 240 MeV/nucleon) one- 
and two-neutron removal reactions on a carbon target. Measurements were made of the inclusive 
cross sections and momentum distributions for the charged residues. Narrow momentum distri- 
butions were observed for one-neutron removal from ^^C and ^"C and two-neutron removal from 
^^C. Two-neutron removal from '^"C resulted in a relatively broad momentum distribution. The 
results are compared with eikonal-model calculations combined with shell-model structure informa- 
tion. The neutron-removal cross sections and associated momentum distributions are calculated for 
transitions to both the particle-bound and particle-unbound final states. The calculations take into 
account the population of the mass A—1 reaction residues, ^~^C, and, following one-neutron emis- 
sion after one-neutron removal, the mass A — 2 two-neutron removal residues, ^~^C. The smaller 
contributions of direct two-neutron removal, that populate the ^~^C residues in a single step, are 
also computed. The data and calculations are shown to be in good overall agreement and consis- 
tent with the predicted shell-model ground state configurations and the one-neutron overlaps with 
low-lying states in ^*~^^C. These suggest significant i^s^/j valence neutron configurations in both 
^"C and ^^C. The results for ^^C strongly support the picture of ^^C as a two-neutron halo nucleus 
with a dominant i^Sj^j ground state configuration. 



I. INTRODUCTION 

Residue momentum distributions following the disso- 
ciation of halo nuclei have long been recognized as probes 
of the spatially extended valence neutron wave functions 
P, Q • Over the last decade, high-energy nucleon removal 
reactions have become a powerful spectroscopic tool to 
explore the structure of nuclei far from stability. Specif- 
ically, the momentum distributions and associated cross 
sections offer a means to probe the active single-particle 
orbitals near the Fermi surface, whereby the shapes of 
the momentum distributions reflect the orbital angular 
momentum of the removed nucleon(s) and the cross sec- 
tions the spectroscopic strength 

Two-nucleon (2iV) removal reactions are challenging 
experimentally and remain less well studied jsj-fllj. Re- 
cently, it has been clarified that 2N removal cross sec- 
tions and the associated residual nucleus momentum dis- 
tributions, following direct removal of two well-bound like 
nucleons, can provide a sensitive probe of rapid struc- 
tural changes |l2l . [l3l |. of spins of final-states [3, [l5l |. 
and also of aspects of two-nucleon spatial correlations 
near the nuclear surface [ll|, [l^- However, a complica- 
tion arises when discussing reactions that result in the 
removal of two weakly-bound (valence) neutrons from a 
neutron-rich, mass A, projectile. The mass A — 2 reac- 



tion residues can now result from two distinct mecha- 
nisms: single-step direct knockout of a pair of neutrons 
and one-neutron knockout followed by neutron emission 
from excited particle-unbound states of the intermediate 
mass A — 1 residue. The latter will be populated with 
cross sections that are typically an order of magnitude 
or more larger than those for the direct pair removal. 
Both processes are discussed quantitatively in the anal- 
yses presented here. 

This article reports the first measurements of one and 
two-neutron removal from i9'20,22(-i secondary beams at 
about 240 MeV/nucleon. The carbon isotopic chain is of 
considerable interest from a structural point of view, as 
it exhibits large odd-even staggering in the one-neutron 
separation energies (see e.g.. Fig. 2 of and weak 

s-wave valence neutron binding and halo formation for 
both [l|. Sin = 1.22(1) MeV, and Sm 
— 0.58(9) MeV. The next heavier odd-A halo candidate, 
^^C, is known to be particle unbound while ^^C is bound 
by S2n = 0.42(94) MeV, which is estimated on the basis 
of systematics [2d|. As such, ^^C may be pictured as a 
three-body {^^C + n + n) Borromean system - having 
no bound two-body subsystems. The structures of the 
i9-22(-^ ground and excited states are considered here by 
the use of the fast neutron removal reaction methodology. 

A recent measurement of the reaction cross section of 



^^C on a proton target [2lj has also drawn much at- 
tention. Its observation of an enhanced cross section 
suggests that ^^C has (a) an extended matter density 
and (b) a two-neutron halo ground state structure, dom- 
inated by a weakly-bound s-wave two-neutron configura- 
tion. Such measurements, carried out over a wide energy 
range, have been successful in deducing effective nuclear 
matter radii and their evolution along isotopic chains. 
However, they cannot provide specific information about 
emerging changes in the microscopic nuclear structures. 

Fast one- and two-neutron removal reactions from the 
lighter neutron-rich carbon isotopes, up to ^^C, have been 
studied experimentally using beryllium and carbon tar- 
gets (see e.g., [isl, [23 - l29| ). Results for neutron removal 
from 18,19q proton target were reported and ana- 

lyzed by Kondo et al. [s^l- The incident beam energies 
of these studies ranged from 45—103 MeV/nucleon. In 
several cases the parallel or transverse momentum dis- 
tributions of the reaction residues were also measured. 
More recently, Ozawa et al. jsH measured one- and two- 
neutron removal reactions of i^'^oc using a proton target 
at 40 A MeV. However, it was difficult to discuss the re- 
actions for due to the low statictics. A systematic 
analysis of many of these results, using the shell model 
and reaction formalism exploited here, was presented in 
Ref. [13. 

As mentioned above, the measurements reported here 
were made at about 240 MeV/nucleon. At such an en- 
ergy domain the underlying assumptions of the sudden 
and eikonal approximations for the reactions are well 
founded. This paper thus presents a quantitative analy- 
sis of the results of both (a) direct one- and two-neutron 
removal, and (b) indirect two-neutron removal. In doing 
so we are able to elucidate the dominant single-particle 
structure of i9'20,22(-i ^^-^^ ^.j^g extent to which these are 
tracked by shell-model calculations at the boundary of 
the p-shell and the lower part of the sd-shell. Only the 
shell- model currently provides this level of information on 
the nucleon single-particle structures consistently along 
an isotopic chain in a form that can be readily used in 
reaction calculations. Thus we employ a consistent shell- 
model and reaction theory description for all of the iso- 
topes considered. 

In Section |ll] we discuss the experimental techniques. 
The results are described in Section IIIII The theoreti- 
cal approach used is outlined in Section |IVl There both 
the reaction theory and the shell-model spectroscopic 
strengths and description of the neutron bound states 
needed for the calculations are discussed. The experi- 
mental results and calculations are compared and dis- 
cussed in Section |Vl The paper then concludes with a 
brief summary. 



II. EXPERIMENTAL SETUP AND DETAILS 

The experiments were performed at the RI Beam fac- 
tory (RIBF) operated by the RIKEN Nishina Center and 
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BigRIPS Secondary target 



FIG. 1: A schematic view of the layout of BigRIPS (F0-F8) 
and the zero degree spectrometer, ZDS (F8-F11), where Fl 
through Fll represent each focal plane. The secondary car- 
bon target was mounted at the achromatic focal plane F8. 
The I8'i9'20c residues were identified using the ZDS. 



the Center for Nuclear Study, University of Tokyo. Sec- 
ondary beams of i9.20,22q -y^rere produced by fragmenta- 
tion of a *^Ga, beam at 345 MeV/nucleon on a 2G-mm- 
thick rotating Be target. The typical "^^Ca primary beam 
intensities (/i) for each setting are listed in Table HI The 
secondary beams were separated using the superconduct- 
ing separator BigRIPS [3.2, .33j whose layout is shown in 
Fig. [T] An achromatic aluminum energy degrader of 15- 
mm-median thickness was located at the dispersive focal 
plane Fl. For the setting to produce the ^^C beam, an 
aluminum degrader of 8-mm-median thickness was also 
installed at the dispersive focal plane F5. The secondary 
beam particles were identified event-by-event by combin- 
ing the measured time-of-flight (TOF) , energy loss {AE) , 
and magnetic rigidity {Bp). The TOF was recorded be- 
tween two plastic scintillators at the achromatic focal 
planes F3 and F7, AE was measured using a plastic 
scintillator at F7, and Bp was determined from a posi- 
tion measurement using a plastic scintillator read out on 
both sides by photomultiplier tubes at F5. An example of 
the particle-identification spectrum is shown in Fig. HJa) 
in terms of the atomic number (Z) and mass-to-charge 
ratio {A/Z), which demonstrates the clear separation of 
^^C. The secondary beam intensities {I2) and the mo- 
mentum spread (AP/P) of the i9^20,22q i^gg^j-Qg ^j-g listed 
in Table m A secondary carbon target (4.02 g/cm^) was 
mounted at the achromatic focal plane F8. The mid- 
target energies (E/A) of i9.20,22(-; ^^^.^ 343, 241, and 
240 MeV/nucleon, respectively, as listed in Table HIl 

The i8,19,20q residues following the reaction were col- 
lected by tuning the rigidity of the zero degree spectrom- 
eter (ZDS) to center the momentum distribution. The 
residues as well as the secondary beam particles were 
identified event-by-event by combining TOF, AE, and 
Bp, where TOF was measured between two plastic scin- 
tillators at the achromatic focal planes F7 and Fll and 
AE was measured using an ionization chamber at Fll. 
The Bp was determined from a position measurement us- 
ing PPACs at the dispersive focal plane F9. The result- 
ing particle- identification spectrum for the '^'^G residues 
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FIG. 2: (color online). Particle-identification spectrum of (a) 
the ^^C beam by BigRIPS, and (b) the residues in the ZDS 
after selecting ^^C ions before the secondary target. 



residues owing to reactions in the target. To do so, total 
reaction cross sections of 1375 mb for ^^C, 1090 mb for 
20c, 1139 mb for ^^C, and 1023 mb for ^^C, which were 
estimated using eikonal calculations, were employed. It 
should be noted that a 10% deviation in these reaction 
cross section values translates into a (essentially negligi- 
ble) 1% deviation of the deduced one- and two- neutron 
removal cross sections. To obtain the residue momentum 
in the projectile rest frame, the Lorentz transformation 
was applied to the laboratory system residue momen- 
tum extracted from the TOF between F8 and Fll. The 
velocity of the projectile was obtained from the TOF be- 
tween F3 and F7. The cross sections c^J'n and the par- 
allel momentum P|| distributions were corrected for the 
transmission efficiency through the ZDS. These were es- 
timated using a Monte Carlo simulation, an angular ac- 
ceptance of the ZDS obtained by calibration runs using 
the secondary beams, and the design value of the momen- 
tum acceptance (AP/P < ±3%) of the ZDS. To obtain 
a higher transmission the momentum acceptance of Bi- 
gRIPS was restricted to AP/P < ±2% in the analysis. 
The transmission was about 90% in all cases. 



is shown in Fig. [H^b). 

The present study measured one- and two-neutron re- 
moval cross sections (o'™^„, x = 1 or 2) and the residue 
parallel momentum (P||) distributions for the channels 
listed in Table HIl The cross sections were derived from 
the numbers of the beam particles counted before the 
secondary target and those of the residues registered at 
the final focal plane Fll of the ZDS, using 



N' 



(1) 



Here Ni {No) represents the number of projectiles and N[ 
{N'^) the number of residues for target-in [i) (target-out 
(o)) runs. Nt and cr/j (cr^) are the number of target nuclei 
per unit area and the reaction cross section of the pro- 
jectile (residue) . The background events were subtracted 
by using target-out runs. Owing to the substantial thick- 
ness of the carbon target, correction had to be applied 
to account for the losses in the number of projectiles and 



III. EXPERIMENTAL RESULTS 

The one- and two-neutron removal cross sections ex- 
tracted here are summarized in Table [IT] together with 
the corresponding mid-target energy of the projectile. 
The cross section for C(^^C, ^®C) is smaller than that 
measured at lower energies. Near 60 MeV/nucleon for 
instance [28|, ct_i„ — 226(65) mb. In part this reflects 
the smaller intrinsic nucleon-nucleon cross sections at 
the energy of the current experiment (being closer to 
the minimum near 300 MeV) as well as changes in the 
real parts of the optical potentials and a reduced diffrac- 
tive breakup contribution for energies in excess of 100 
MeV/nucleon For ^Oc, the yield of ^^C residues is 
smaller than that for ^^C, reflecting the available energy 
windows for final states below the first neutron thresholds 
in i^c (S'„ = 4.18 MeV) and ^^C (5'„ = 0.58(9) MeV). 
For 22 C we note an enhanced cross section compared to 
2^0 which, as we will show, reflects the two-neutron halo 
character of 22c. 



TABLE I: The typical ''**Ca primary beam intensity (/i) for 
each setting and typical secondary beam intensity (/2). Also 
tabulated are the momentum spread ( AP/ P) of the secondary 
beams. 

Secondary beam I\ (pnA) I2 (particles/s) AP/P 



19^ 
20c 
22 ri 



^4 
^80 



X 10'' 
^6 X 10^ 
^6 



±2% 
±3% 
±3% 



TABLE II: The one- and two-neutron removal cross sections 
for each reaction channel at the mid-target energies (E/A) 
shown. The widths (FWHM) of the momentum distributions, 
after unfolding the experimental resolutions, are also shown. 



Reaction 



E/A (MeV) 



(mb) FWHM (MeV/e) 



C("C. 

cC°c. 



i«C) 
«C) 
i«C) 
2°C) 



243 
241 
241 
240 



163(12) 
58(5) 
155(25) 
266(19) 



73(5) 
46(13) 
204(11) 

76(8) 
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The measured momentum distributions are shown in 
Fig.El That for (^^C, 2°C) was obtained for the first time 
here. The shapes of the distributions for ("^^C, ^*C), (^°C, 
^^C), and (^^C, ^°C) are well reproduced by a Lorentzian 
distribution convoluted with the experimental resolutions 
of 23, 23, and 27 MeV/c, respectively. For {^°C, ^^C), 
the distribution is better reproduced by a Gaussian con- 
voluted with the experimental resolution of 28 MeV/c. 
The widths (FWHM) after accounting for the experimen- 
tal resolutions are listed in Table HT] 

The distributions for (^^C, ^^C) and (^^C, ^OC) are rel- 
atively narrow, with widths of 73(5) and 76(8) MeV/c. 
Qualitatively, these arise from the weakly-bound one and 
two-neutron halo-like nature of these nuclei and the role 
of s-wave valence neutron(s) in their ground state struc- 
ture. It is interesting that the measured distribution for 
(^"C, ^^C) is also extremely narrow. Qualitatively, this is 
the result of the reaction probing directly the com- 
ponent of the ^°C(0+) ground state wave function, as is 
required to populate the bound "'^^Cg.s.(l/2^) halo-state 
residues. The quantitative analyses for these momentum 
distributions and their cross sections are detailed in the 
following sections. 



IV. REACTION ANALYSIS 
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A. Reaction theory 



FIG. 3: The momentum distributions for (a) C("C, ^^C), 
(b) CC"C, ^»C), (c) C(2°C, ^«C), and (d) CC^C, ^"C). The 
shapes of the distributions for cases (a), (b), and (d) are well 
reproduced by a Lorentzian distribution convoluted with the 
experimental resolution. For the (^°C, ^*C), case (c), the 
distribution is better reproduced by a Gaussian. The distri- 
bution shapes were fitted over the momentum ranges —200 to 
+200 MeV/c. 



We adopt an eikonal model description of the reaction 
mechanisms. Given a nucleon- or nucleus-target inter- 
action description, the eikonal approximation has been 
shown (35} to provide a rather accurate description of 
the elastic 5'-matrix and derived observables for incident 
projectile energies of order 50 MeV/nucleon and greater. 
As noted earlier, at the energy of the current experi- 
ments (about 240 MeV/nucleon) the sudden and eikonal 
approximations of the reaction model are very accurate. 

The removal reaction cross sections for one-neutron 
knockout to a given final state, with spin-parity , are 
calculated using Q 



E 

nlj 



A 



A-l 



C^S{J\nlj)a,^{nt3,S\ 



'Off\ 



(2) 



where the C^S are the shell model spectroscopic factors 
and the single-particle cross section (Jsp is calculated us- 
ing the eikonal model assuming unit spectroscopic fac- 
tor. The quantum numbers of the removed neutron are 
denoted by n€j and 5*°* is the effective separation energy 
for the removal of the neutron leaving the residue in the 
given final state. The single-particle cross sections, Csp, 
include the contributions from both the stripping (inelas- 
tic breakup) and diffractive dissociation (elastic breakup) 
mechanisms. Details of calculations of these two distinct 
and (incoherent) additive contributions can be found in 
Ref. [3I. 



In direct two-neutron removal the theoretical cross 
sections do not in general (e.g., when there are sev- 
eral active orbitals) factorize into a structural (spectro- 
scopic) factor and a single-particle cross section. The 
cross sections involve coherent contributions from all ac- 
tive shell-model two-nucleon configurations with non- 
vanishing two- nucleon amplitudes (TNA) . Details of their 
definitions and the phase conventions used can be found 
in Ref. [9]. Here we will calculate the single-step direct 
two-neutron removal yields arising from both (a) two- 
neutron stripping and (b) one neutron being stripped and 
the second being elastically dissociated (diffracted) fll| . 
Since these direct two-neutron removal cross sections are 
small compared to the cross sections arising from single- 
neutron removal, we do not expand upon these formal 
aspects here. Full details of the necessary eikonal formal- 
ism, as is applied to direct two-nucleon removal events, 
can be found in Refs. [ol. [T]|. 

For both the one- and two-neutron removal calcula- 
tions, the required neutron- and residue-target elastic S- 
matrices were calculated using the static density limit 
of the eikonal model, e.g., |37| . also known as the tpp 
limit of the Glauber multiple scattering series. That 
is, we used the single-folding model (nucleon-target) 
and double-folding model (residue-target) for the ab- 
sorptive optical model interactions with the carbon tar- 
get. The inputs needed were the residue and target 
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point neutron and proton one-body densities and an ef- 
fective nucleon-nucleon (NN) interaction. The densities 
of the i8.i9,20,2iQ j-gsijues were estimated from spherical 
Skyrme Hartree-Fock (HF) calculations using the SkX 
interaction [sl]. The HF calculates the (experimentally 
unbound) ^^C case to be weakly bound. We used this 
bound density for calculations of the ^-'^C-target optical 
potential (and its ^-matrix) in the localized region where 
they make grazing collisions. 

All calculations assumed the following. The density 
of the carbon target nuclei was taken to be of Gaus- 
sian form with a point-nucleon root-mean squared radius 
of 2.32 fm. A zero-range NN effective interaction was 
used with its strength calculated from the free neutron- 
neutron and neutron-proton cross sections at the beam 
energy and from the real-to-imaginary ratios of the NN 
forward scattering amplitudes at 240 MeV, interpolated 
(using a polynomial fit) from the values tabulated by Ray 
|39| . The use of these inputs, as a standard parameter set 
in the eikonal reaction model, was shown to accurately 
reproduce the recently-measured ^oj ratios of the diffrac- 
tion to stripping reaction mechanism yields in the cases 
of *B and induced proton-removal reactions. A recent 
careful analysis by Bertulani and De Conti [4l| confirms 
that corrections to this adopted procedure, due to Pauli 
blocking corrections to the NN effective interaction, are 
negligible at the energy of the present study. 

We assume here that the heavy residue-target inter- 
actions and their S'-matrices are diagonal with respect 
to the different final states of the residue, and thus that 
there is no reaction-induced dynamical excitation of the 
residues during the collision. For the odd-A carbon 
projectiles, where different neutron orbitals (nij) may 
then contribute to a given J'^ final state, this implies 
that these different n£j contributions are incoherent and 
should be summed. 



using a standard Woods-Saxon potential geometry with a 
reduced radius parameter tq — 1.25 fm and a diffuseness 
of flQ — 0.7 fm. 

The depths of the potential wells in each instance were 
adjusted to reproduce the neutron separation energy, tak- 
ing into account the excitation energy of the final state. 
The ground state to ground state one- and two-neutron 
separation energies were taken from the 2003 mass evalu- 
ation [l^l • The shell-model excitation energies of Tables 
mil and HVl were then used (without any adjustment) for 
all non-ground-state transition calculations. For the di- 
rect two-neutron removal calculations each neutron was 
assumed to have a separation energy of half that for the 
two neutrons to the final state of interest. 



V. ANALYSIS AND DISCUSSION 

We discuss the results for the inclusive cross sections 
for one- and two- neutron removal from the 19.20,22 q jgQ_ 
topes. We consider in detail the calculated contributions 
from both indirect and direct two-nucleon removal. The 
measured and calculated inclusive parallel momentum 
distributions are also discussed. In all cases these are 
shown in the projectile rest frame. In the comparisons 
with the data the theoretical momentum distributions, 
calculated using the stripping mechanism, have been con- 
voluted with the Gaussian experimental resolution given 
in Section IIIII and then scaled to the measured inclusive 
cross sections. Further discussion of the calculations of 
the parallel momentum distributions in the cases of the 
transitions to unbound final states is included when dis- 
cussing these; i.e., for the 20,22 q cases. 

A. One-neutron removal reactions 



B. Shell-model calculations and overlaps 

A consistent set of shell-model calculations were used 
for the required level energies, spectroscopic factors and 
two-nucleon amplitudes. These were performed using the 
code OXBASH 42]. The calculations used the Warbur- 
ton and Brown (WBP) effective interaction [4^ in a psd- 
model space truncated to allow Ohu and Ihtj excitations. 
The small center of mass correction factor, [A/(A — 1)]^, 
shown in Eq. where TV is the principal oscillator 
model quantum number of the orbit of the removed nu- 
cleon ;44J, was applied to the (fixed center) shell- model 
spectroscopic factors of all of the single-neutron removal 
calculations. 

The low-lying shell-model levels and spectroscopic fac- 
tors for the reaction products, i8,i9,20,21q^ g^^^ collected 
in Tables IIIII and IIVI The wave function of each of 
the removed-neutron bound states was calculated in a 
Woods-Saxon potential well of a fixed geometry. Follow- 
ing Ref . 17[ : these radial wave functions were calculated 



We first discuss the individual and the inclusive one- 
neutron removal cross sections to bound, and un- 
bound (neutron emitting), <J^\^(^^y states of the mass 
A — 1 isotopes. The experimental and theoretical results 
of the present work are collected in Tables IIIII and IIVI 
Also tabulated are the details of the relevant shell-model 
states, their energies, spins, parities and spectroscopic 
factors, C^S. The overall ratios of the measured to the 
calculated inclusive one-neutron removal cross sections, 
Rs = (j'^in/cr^-in: a-i'e also shown in the tables for each 
bound and unbound final states data set. 



1. Results for C 

The case of provides a valuable link to the earlier 
work at lower energies, summarized in [Tt} , and the re- 
lated and more exclusive results using neutron knockout 
from a proton target [30| . From the present work, calcu- 
lated exclusive and experimental inclusive one-neutron 
removal yields from the -'^^C(l/2+) ground state, with 
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ground-state separation energy S'i„(-^^C) = 0.58 MeV, 
are shown in Table IIIII The theoretical cross sections 
are shown for the six predicted positive parity ^^C final 
states. 

In the case of (^^C, ^^C) the WBP shell model calcu- 
lation places several final states near to or between the 
one- and two- neutron threshold energies of 4.18 MeV and 
4.91 MeV, respectively, for ^^C. Specifically, the third 
2;^ and first 3^ states at 4.915 and 4.975 MeV have 
significant spectroscopic strengths and associated cross 
sections. Experimentally, recent work of Kondo et al. 
[30l |. on neutron knockout from ^^C on a proton target, 
observed gamma-rays from a (2+,3"'") excited state (or 
states) near 4.0 MeV, the associated "'^^C transverse mo- 
mentum distribution being characteristic of an £ = 2 neu- 
tron removal. The earlier (^C, ^*C) inclusive data anal- 
yses of Maddalena et al. [2J| and Simpson and Tostevin 
[17| also assumed these 2;^ and 3^ states near 4.9 MeV 
were neutron bound, citing the results of shell-model cal- 
culations using a modified version of the WBT interac- 
tion [H. 

We have also calculated the inclusive parallel momen- 
tum distributions, to bound final states, as the sum of 
the distributions to these individual final states weighted 
by the (J^\n shown in Table Hill Fig. 2] shows the experi- 
mental (^^C, ^^C) inclusive parallel momentum distribu- 
tion and also those calculated. In all cases the theoreti- 
cal momentum distribution curves are normalised to the 
measured inclusive cross section. We show the results 
obtained by (a) assuming that the 2^ and 3^ states are 
unbound (dashed curve), having an inclusive cross sec- 
tion of 100.2 mb, and (b) assuming that the 2^ and 3^ 
states are bound (solid curve), resulting in an inclusive 
cross section of 164.4 mb. The experimental cross sec- 
tion, in Table Hill is 163(12) mb. The comparison with 
the present momentum distribution data, in particular, 
provides us with rather compelling evidence for the hy- 
pothesis (b), that the 2;^ and 3]^ states are bound. 



2. Results for C 

The predicted ^^C shell-model final states and the cal- 
culated and experimental one-neutron removal cross sec- 
tions from the ^°C(0^) ground state, with separation en- 
ergy 5i„(20C) = 2.90 MeV are collected in Tables Imj 
and IIVI There is only very incomplete experimental 
information on the low-lying excited state spectrum of 
^^C. Using Coulomb dissociation the ^^C ground state 
has been unambiguously identified as a 1/2+ s-wave halo 
state with weak binding [l^. The evaluated ^^C first 
neutron threshold is at 0.58(9) MeV ^2d\. An unbound 
excited 5/2+ state with = 1.46(10) MeV has also been 
clearly identified (46| using inverse-kinematics proton in- 
elastic scattering from -^^C. Stanoiu et al. reported a 
201(15) keV gamma-ray transition in ^^C with in-beam 
7-ray spectroscopy. Using inverse-kinematics proton in- 
elastic scattering Elekes et al. 47[ also reported evidence 
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FIG. 4: Measured inclusive parallel momentum distribution 
of i*C, following one-neutron removal from ^^C on a carbon 
target at 243 MeV/nucleon compared to the theoretical calcu- 
lations. The solid line includes contributions from the 2^ and 
3+ shell-model states of ^*C, assumed bound; see also Table 
IIIII The dashed line shows the results when assuming that 
the 2 J and 3+ states are unbound. Here, and in Fig. [S]|S1 the 
theoretical distributions have been convoluted with the exper- 
imental resolution and normalized to the measured inclusive 
cross section. 



of two gamma-ray transitions, with energies 72(4) and 
197(6) keV, from two bound ^^C excited states. While in 
both of these cases the transition energy (near 200 keV) 
is close to that of a predicted 5/2+ shell-model bound 
excited state, we will show that the present experimental 
data and analysis do not support such a strong transition 
to a 5/2+ bound ^^C excited state. 

Table Hill shows the cross section for the 1/2+ shell- 
model ground state transition. The measured cross sec- 
tion, of 58(5) mb, and parallel momentum distribution 
to a bound ^^C final state, Fig.[5l are consistent with the 
theoretical expectations for the removal of a 2si/2 neu- 
tron (solid curve) with the tabulated shell-model spec- 
troscopic factor of near to unity. It is likely that the 
1/2+ ground state of ^^C is the only bound state of this 
system. If any bound excited state exists, the only possi- 
bility seems to be the first 3/2+ state predicted at 0.624 
MeV, which would add the cross section of only 7.20 mb. 

Although the bound 5/2+ state is unplausible due to 
the observed small cross section, we attempt to estimate 
an upper limit on possible bound d-state strength below. 
Table |TV| shows the results for the cross sections lead- 
ing to the excited ^^C shell-model final states. The shell 
model predicts seven such excited states with significant 
spectroscopic factors below the ^^C neutron threshold of 
4.18 MeV. Given these cross sections we note, from Table 
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FIG. 5: Measured inclusive parallel momentum distribution 
of i^C, following one-neutron removal from on a car- 
bon target at 241 MeV/nucleon compared to the theoreti- 
cal calculations. The solid curve assumes that only the 1/2+ 
shell-model ground state transition (2si/2 neutron removal) 
is bound. The long-dashed, short-dashed, and dot-dashed 
curves result if one assumes that 0.5, 1.0 or 1.5 units of 
spectroscopic strength also leads to bound final states. 



FIG. 6: Measured inclusive parallel momentum distribution 
of ^*C, following two-neutron removal from ^"C on a carbon 
target at 241 MeV/nucleon compared to the theoretical cal- 
culations. The theoretical curves are the weighted sum of the 
exclusive calculations of the unbound ^^C states, see text. Re- 
coil effects associated with the neutron emission are included 
assuming the most important contributions come from states 
with e* of 1.0 MeV (dashed curve) and 2.0 MeV (solid curve). 



IIV[ that one unit of the first excited state 1^5/2 spectro- 
scopic strength makes a contribution of 30.5 mb to the 
theoretical cross section. The calculated 1/2^ ground 
state cross section is 58.92 mb. Thus, if there was also 
1^5/2 strength to bound state(s), with a summed spec- 
troscopic strength of 0.5, 1.0, or 1.5 units, the theoretical 
cross section to bound final states would increase to 74.2, 
89.4, or 104.7 mb, respectively, well in excess of the mea- 
sured value of 58(5) mb. The corresponding effects of 
such bound lc?5/2 strength on the shapes of the calcu- 
lated parallel momentum distributions are shown in 
Fig. [5] by the long-dashed, short-dashed, and dot-dashed 
curves. Here, each curve is normalized to the experi- 
mental cross section of 58(5) mb. We conclude from this 
comparison that the majority of the strength that the 
shell- model attributes to the 190 keV lc?5/2 state is in 
fact unbound. Based on Fig. [5] and the measured cross 
section to bound ^^C, we estimate that 0.5 units or less 
of bound 1(^5/2 strength might be accommodated by the 
present data set. 

Our assumption, in Table ITVl is that all of the excited 
"'^^C shell-model states are unbound and that these un- 
bound states will decay by neutron emission to ^*C. In 
this and the following case of ^^C these unbound mass 
A—1 excited state cross sections are large. For such un- 
bound final state cases our one-neutron removal model 
calculates the exclusive parallel momentum distributions 



of the (weakly) unbound ^^C and ^^C residues in the 
original projectile rest frame. The subsequent in-flight 
neutron emission from these excited states will generate 
additional (recoil) broadening of the momentum distri- 
butions of the observed mass A — 2 residues, i.e., and 
^°C. The degree of broadening will be dependent on the 
continuum energy of the unbound, mass A — I interme- 
diate state, denoted e*. 

We estimate the effect of this recoil. We assume that, 
in the rest frame of the unbound, mass A — 1 state, with 
its given continuum energy e* , the mass A — 2 residue (in 
its ground state) and the decay neutron separate isotropi- 
cally. The momentum p of the heavy decay residue in this 
frame satisfies = 2^e* , where ji is the A — 2 residue- 
neutron reduced mass. The assumption that this two- 
body decay is isotropic then requires that the calculated 
parallel momentum distributions of the unbound mass 
A — 1 fragments must be convoluted with a rectangu- 
lar distribution, of unit integral and total width 2p, to 
derive the mass A — 2 fragment parallel momentum dis- 
tributions. This is done for the theoretical distributions 
shown in Figs. [6] and[71 

Table IIVI shows that these indirect two-neutron re- 
moval cross sections arise predominantly from interme- 
diate states of small e* . For the case, a cross section 
of 160 mb is predicted to arise from the first three shell- 
model excited states having < 1 MeV in -"^^C. How- 
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ever, as was discussed above, these shell model energies 
are not sufhciently accurate. The WBP interaction pre- 
dicts the first excited 5/2^ state to be bound with a large 
spectroscopic factor, whereas the calculated (y^\n to this 
state and the measured cross section and momentum dis- 
tribution to the ^^C ground state exclude this possibility. 
Experimentally a ^^C excited 5/2+ state has been clearly 
identified at Ey, = 1.46(10) MeV by Satou et al. 46] in 
proton inelastic scattering from ^^C, i.e., with e* = 0.88 
MeV. The present inclusive data do not permit a more 
detailed analysis of this excited state or of the predicted 
shell-model strength distributions. Our analysis shows, 
however, that the present data are consistent with an in- 
tegrated 5/2+ strength of about 4 units leading to the 
"'^^C continuum, as is given by the shell model. 

In the absence of more complete information, and to 
assess this recoil sensitivity, we calculate the evaporation 
recoil effects assuming that the most important contri- 
butions arise from ^^C intermediate states with (i) e* — 
1.0 MeV and (ii) e* = 2.0 MeV. The results are shown 
in Fig. [5] by the dashed line and solid line, respectively. 
Both outcomes are consistent with the experimental ^^C 
residue momentum distribution. We conclude from this 
agreement of the inclusive cross section and the momen- 
tum distribution that the 5/2+ shell- model excited state 
is very probably unbound. Since, in this case, the shell 
model appears to systematically produce states with too 
small an excitation energy, the effective neutron separa- 
tion energies will also be underestimated and, in turn, the 
theoretical removal cross sections slightly overestimated. 
This is also suggested by the smaller Rs — ^^-2n/ '^'^\n(e) 
deduced for this isotope. We do not attempt to make any 
parameter adjustments to compensate for this (small) ef- 
fect. 



3. Results for C 

Here all final states of the ^^C one- neutron removal 
residues are particle unbound. The calculated exclusive 
(and inclusive) and experimental inclusive one-neutron 
removal yields are collected in the Table IIVI for the pre- 
dicted shell-model states of ^^C, that decay by neutron 
emission to ^"^C. 

Very little is known about these isotopes. Both the 
one- and two- neutron separation energies from ^^C are 
only poorly determined and so we are guided by the 
2003 mass evaluation That is S2n{^'^G) = 0.42(94) 
MeV and Sin{^^C) = -0.33(56) MeV, with large uncer- 
tainties. Thus, the ground state of ^^C was assumed to 
be produced at a continuum energy of e* = 0.30 MeV 
after neutron removal with ground state separation en- 
ergy S'i„(^^C) = 0.70 MeV. As was discussed for the 
projectile case, the inclusive (unbound) ^^C momentum 
distribution is calculated as the weighted sum of the mo- 
mentum distributions to the individual final states with 
the cr^\n(e) shown in Table ITVl The neutron emission re- 
coil broadening is included for each final state according 
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FIG. 7: Measured inclusive parallel momentum distribution 
of ■^"C, following two-neutron removal from ^^C on a carbon 
target at 240 MeV/nucleon compared to the theoretical calcu- 
lations. The solid curve is the weighted sum of the exclusive 
calculations for the unbound ^^C states, see text. The dashed 
and dot-dashed curves show the contributions from knockout 
via the 1/2+ and 5/2+ unbound '^^C intermediate states, re- 
spectively. The recoil broadening arising from neutron emis- 
sion from these unbound intermediate states is folded in. 



to its e* value, i.e., e* = E^ + Q.'iO MeV, prior to this sum 
being performed. Three final states are predicted below 
the first neutron threshold of 2.90 MeV. 

Table lTVl shows that, based on the shell-model, the first 
two final states each contribute almost half of the inclu- 
sive one-neutron removal cross section. These states are a 
1/2+ ground state, with spectroscopic factor C'^S = 1.4, 
and a 5/2+ neutron- hole state at = 1.11 MeV, with 
C'^S = 4.2. The associated measured and theoretical in- 
clusive parallel momentum distributions (convoluted 
with the experimental resolution of 27 MeV/c) arc com- 
pared in Fig. [7l The individual contributions from the 
two dominant shell-model final states are also shown in 
the figure. The agreement with the data is very good, 
providing strong support for the weakly-bound, vs\i2 

character for the ^^C ground state. This result is con- 
sistent with the recent interaction cross section measure- 
ment and associated analysis of Ref . [2l[ , that is sugges- 
tive of an extended ^^C matter density. 

Currently, the ^^C two-neutron separation energy, 
S'2„(22C) = 0.42(94) MeV, has a significant uncertainty. 
Hence, we consider the sensitivity of the theoretical in- 
clusive ^''C production cross section and momentum dis- 
tribution to the value assumed. Fig. [5] shows the cal- 
culated momentum distributions when assuming a ^^C 
two-neutron separation energy of S2n(^'^G) = 0.40 MeV 
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FIG. 8: Measured inclusive parallel momentum distribution 
of ^"C, following two-neutron removal from ^^C on a carbon 
target at 240 MeV/nucleon compared to the theoretical cal- 
culations. The theoretical curves are the inclusive cross sec- 
tions calculated assuming ^^C two-neutron separation ener- 
gies S2nC^C) = 0.40 MeV (solid curve), 0.70 MeV (dashed 
curve), and 1.20 MeV (dot-dashed curve). The curves include 
the recoil broadening arising from the neutron decay of the 
unbound ^^C intermediate states. 

(solid curve), 0.70 MeV (dashed curve), and 1.20 MeV 
(dot-dashed curve). In these calculations we continue to 
assume that the ground state of ^^C is at a continuum en- 
ergy of 0.30 MeV, hence the ground state to ground state 
one-neutron separation energy of ^^C is S'i„(^^C) = 0.70, 
1.00, and 1.50 MeV in these cases. The curves have been 
convoluted with the experimental resolution of 27 MeV/c 
and also include the recoil broadening arising from the 
neutron decay of the unbound ^^C states ( Table ITV)) . The 
increasing separation energies reduce the corresponding 
inclusive cross sections: 283, 257, and 227 mb for the 
52„(2^C) = 0.40, 0.70, and 1.20 MeV, respectively We 
note that, owing to the relative insensitivity of our calcu- 
lated cross sections and momentum distributions to the 
current experimental uncertainty in S'2n(^^CI), the data 
of the present work do not determine or place a signifi- 
cant constraint upon this value. We are able to conclude, 
however, that the sensitivity to the underlying structure 
of ^^C, specifically of the approximately equal contribu- 
tions of the 1/2^ and 5/2'^ transitions to the measured 
inclusive cross section, is robust. 

B. Direct two-neutron removal 

We summarize only briefly the calculated exclusive 
and inclusive direct two-neutron removal cross sections. 



cr_2n{d)' from and ^^C to bound states of the mass 
A — 2 isotopes. These results are collected in Table |Vl 
computed based on the WBP shell-model two nucleon 
amplitudes, TNA. As noted earlier, as these calculated 
cross sections were both expected and found to be small, 
we will not enter into an extended discussion and details 
of the calculations. The descriptions of the nucleon over- 
lap functions used and the construction of the residue- 
and neutron-target S-matrices are the same as for the 
one-neutron removal analysis. For the full details and 
the formalism of the exclusive cross sections (and their 
momentum distributions) the reader is referred to recent 
references [TTI. [T5j|. 

For projectiles six states below the neutron thresh- 
old in ^^C have appreciable TNA. These states include 
the 2^ and 3^ states proposed as being bound from the 
one- neutron removal analysis (Section lV A ip . The inclu- 
sive direct two-neutron removal cross section is calculated 
to be 14.6 mb. For ^^C projectiles, just two states below 
the neutron threshold in have appreciable TNA and 
the direct two-neutron removal cross section is now 12.1 
mb. These numbers are to be compared with those for 
the indirect two- neutron removal paths that predict cross 
sections of 191.2 and 283.0 mb, respectively. In addition 
we note that, in the case of the removal of strongly-bound 
two-neutron pairs, these calculated direct two-nucleon re- 
moval cross sections typically overestimate the measured 
cross sections with R,{2N)''= cr''_^i^/a%^^^^ « 0.5 [HI- 
Thus, as was found in the earlier study of the lighter car- 
bon isotopes p/7|], the direct pathways enter at about an 
8% level. Since we are unable to distinguish these direct 
events with the current experimental setup, they cannot 
be elucidated or exploited further here. 



VI. CONCLUSIONS 

In summary, the structure of the most neutron-rich 
carbon isotopes i9.20,22(-< j^g^g been investigated using 
high-energy (about 240 MeV/nucleon) single and two- 
neutron removal reactions on a carbon target. Narrow 
momentum distributions were observed for one-neutron 
removal from ^^C and and two-neutron removal from 
^^C. A much broader momentum distribution was found 
in the case of two-neutron removal from ^'^C. 

The measured cross sections and momentum distri- 
butions were interpreted in the light of eikonal reac- 
tion model calculations for single-neutron knockout com- 
bined with structural input derived from psd shell-model 
calculations employing the WBP interaction. The two- 
neutron removal cross sections were calculated by consid- 
ering (a) the removal of one neutron to unbound states 
in the A — 1 daughter, with the assumption that these 
unbound intermediate states decay by neutron emission 
to bound states in the mass A — 2 residue and (b) direct 
two-neutron removal. 

In the case of C(^^C,^°C), the cross section and mo- 
mentum distribution are consistent with the existence 
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of a two-neutron halo with a dominant vs^^^ configu- 
ration in ^^C. The very narrow momentum distribution 
and relatively low cross section for C(^°C/^C) arises as 
the single-neutron removal to the ^^C(l/2+) ground state 
probes specifically the significant vs\^2 component of the 
ground state. The 2^(^5/2 component results in the 
population of unbound states in ^^C that neutron de- 
cay to ^®C. The narrow momentum distribution and en- 
hanced cross section for C(^^C/^C) are consistent with 
the well developed 1^51/2 halo of ^^C. 

Overall, the calculated cross sections agreed well with 
those measured. In particular, in the cases of C(^^C,^^C), 
C(2"C,i9C), and C(22C,2"C), Rs = '^-LZ-m was close to 
unity and consistent with systematics [7[ . Combined with 
the good agreement for the momentum distributions, it 
may be seen that the shell-model has predictive power 
in this region and provides a good overall description of 
level positions and their spectroscopic strengths. 

In the case of weakly-bound two-neutron removal, 
fully-correlated, direct removal cross sections were also 
calculated. It was shown that two-neutron removal is 
dominated by the two-step process; that is, one-neutron 
removal followed by neutron decay of the unbound in- 
termediate state(s). In the cases of two-neutron removal 
from 20,22(-i^ j-j^g direct removal contribution was com- 
puted to be of order 8% or less of the two-step pro- 
cess; consistent with an earlier work [17|. Although this 
makes the identification of these direct two-nucleon re- 



moval events challenging, kinematically complete mea- 
surements of 20,22q breakup are expected to be possible 
in the near future and will help clarify and quantify this 
process. 
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TABLE III: Results for one-neutron removal reactions from ^'"'"'^''C. Tabulated are the one-neutron removal cross sections to 
assumed bound shell-model states near and below the neutron thresholds in the mass A-1 systems, i*'"C, of 4.18 and 0.58 
MeV, respectively (see also the footnotes). The final theoretical cross sections, a-*^i„, include the center-of-mass correction 
factor [A/ {A - 1)]^. The errors shown for the ratio of cross sections, Rs — a'^[^/a^]^i„, reflect only the errors quoted on the 
measurements. 

Reaction E^jMeV) J" £ a,p (mb) C'^S (mb) (mb) Rs 



(i^C(l/2+),i»C(J")) 
5i„(i^C) = 0.58 MeV 



163(12) 0.99(7) 

(2«C(0+),"C(J")) 

SmCC) = 2.90 MeV Inclusive 58.92 58(5) 0.98(8) 



0.000 







104.7 


0.580 


67.63 


2.114 


2^ 


2 


29.91 


0.470 


15.67 


3.639 




2 


25.91 


0.104 


3.00 


3.988 







39.35 


0.319 


13.97 


4.915" 


3^ 


2 


23.60 


1.523 


40.04 


4.975" 


2+ 


2 


23.50 


0.922 


24.15 




Inclusive 








164.5 


O.OOO'' 


1/2+ 





48.37 


1.099 


58.92 




Inclusive 








58.92 



"The 2J and 3+ states at 4.915 and 4.975 MeV are assumed 
to be bound (see Section I V A 111 . 

'There is no evidence from the present work that the i?x = 0.190 
MeV, 5/ 2+ s hell- model state in -"^^C is bound. This state is included 
in Table ITVl and is treated as unbound. 



TABLE IV: Results for the indirect two-neutron removal reaction cross sections. Tabulated are the one-neutron removal cross 
sections to all predicted unbound A — 1-body shell-model states with energies below the neutron threshold of the mass A — 2 



"-ln(e) 
th 

Iri(e) 5 



respectively. The final theoretical cross sections, 
shown on the ratio of cross sections, Rs — fl'2^/0" 

values do not include direct two-neutron breakup events, and the Rs values represent upper limits. 



include the center of mass correction factor [A/{A — 1)]^ 
reflect only the errors quoted on the measurements 



The errors 



The 



-ln(e) 



Reaction 




Sx(McV) 


J" 


£ 


asp (mb) 


C^S 


^-\nM (mb) 




Rs 




,i«C(J^)) 


0.190° 


5/2+ 


2 


27.50 


3.649 


111.17 








= 2.90 MeV 


0.624 


3/2+ 


2 


26.34 


0.247 


7.20 










0.927 


i/2r 


1 


26.46 


1.426 


41.79 










1.541 


5/2+ 


2 


24.31 


0.282 


7.59 










2.417 


3/2r 


1 


22.27 


0.689 


17.00 










3.284 


3/2+ 


2 


21.50 


0.191 


4.56 










3.717 


1/2J 





30.53 


0.055 


1.86 












Inclusive 








191.2 


155(25) 


< 0.81(13) 




''C(J-)) 


0.000 


1/2+ 





89.35 


1.403 


137.55 








= 0.70 MeV 


1.109 


5/2+ 


2 


29.39 


4.212 


135.87 










2.191 


3/2+ 


2 


25.44 


0.342 


9.55 












Inclusive 








283.0 


266(19) 


< 0.94(7) 


"There is no evidence from the present work that 5/2+ 


shell-model 











state in is bound. It is assumed to be unbound. 
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TABLE V: Theoretical results for the direct two-neutron removal reaction cross sections, o-^J^2nid) ■ Tabulated are the two-neutron 



removal cross sections to all predicted shell-model states below the neutron threshold in the mass A~2 systems, ^^'^^C (4.18 
and 2.90 MeV, respectively). 



Reaction 


-C/x(Me V ) 


ITT 

J 


^th /™U\ 


(2''C(0+),i*C(J")) 


0.000 


0^ 


5.66 


S2nC°C) = 3.51 MeV 


2.114 


2f 


4.00 




3.639 


2^ 


0.53 




3.988 




0.36 




4.915" 




1.98 




4.975" 




2.10 






Inclusive 


14.6 


(22C(0+),2«C(J")) 


0.000 


ot 


5.32 


S2n{^^C) = 0.40 MeV 


2.102 




6.81 






Inclusive 


12.1 



"The 2;^ and 3^ st ates at 4.915 and 4.975 MeV are assumed 
to be bound (see Section FV A 11 1. 



